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Abstract — In this paper, we present a simple teciinique to 
approximate the performance union bound of a punctured 
turbo code. The bound approximation exploits only those terms 
of the transfer function that have a major impact on the 
overall performance. We revisit the structure of the constituent 
convolutional encoder and we develop a rapid method to calculate 
the most signiflcant terms of the transfer function of a turbo 
encoder. We demonstrate that, for a large interleaver size, 
this approximation is very accurate. Furthermore, we apply 
our proposed method to a family of punctured turbo codes, 
which we call pseudo-randomly punctured codes. We conclude 
by emphasizing the benefits of our approach compared to 
those employed previously. We also highlight the advantages of 
pseudo-random puncturing over other puncturing schemes. 

I. Introduction 

Turbo codes, originally conceived by Berrou et al. [1] 
are widely known for their astonishing performance on the 
additive white Gaussian noise (AWGN) channel. Methods to 
evaluate an upper bound on the bit error probability (BEP) of 
a parallel-concatenated coding scheme have been proposed by 
Divsalar et al. [2] as well as Benedetto and Montorsi [3]. In 
addition, guidelines for the optimal design of the constituent 
convolutional codes were presented in [4]. 

The rate of a turbo code can be increased by puncturing 
the outputs of the turbo encoder. Guidelines and design 
considerations for punctured turbo codes have been derived 
by analytical [5]-[7] as well as simulation-based approaches 
[8], [9]. Upper bounds on the bit error probability (BEP) can 
be easily evaluated based on the techniques presented in [7] 
and [10]. However, computation of the upper bound can be 
complex and time-consuming, when a large interleaver size 
and certain puncturing patterns are used. 

The motivation for this paper is to derive simple expressions 
for the calculation of the dominant term of the performance 
union bound for punctured parallel concatenated convolutional 
codes (PCCCs). Previously, complex approaches based on the 
full transfer function of each constituent code, have been used. 
In Section we demonstrate that for a large interleaver size, 
the dominant term can be used as an accurate approximation 
of the overall performance union bound. In Section |lll] we 
analyze the properties of constituent convolutional encoders so 
as to obtain exact expressions for the dominant term. A case 
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Study considering pseudo-random puncturing is presented in 
Section |IV] and the paper concludes with a summary of the 
main contributions. 

II. An Upper Bound to the Error Probability of 
Punctured Turbo Codes and its Approximation 

Turbo codes, in the form of rate- 1/3 PCCCs, consist of two 
rate- 1/2 recursive systematic convolutional (RSC) encoders 
separated by an interleaver of size N [1]. The information 
bits are input to the first constituent RSC encoder, while an 
interleaved version of the information bits are input to the 
second RSC encoder. The output of the turbo encoder consists 
of the systematic bits of the first encoder, which are identical to 
the information bits, the parity-check bits of the first encoder 
and the parity-check bits of the second encoder 

Rates higher than 1/3 can be obtained by periodic 
elimination of specific codeword bits from the output of 
a rate- 1/3 turbo encoder Punctured codes are classified as 
systematic (S), partially systematic (PS) or non-systematic 
(NS) depending on whether all, some or none of their 
systematic bits are transmitted [9]. Note that a punctured 
PCCC can also be seen as a PCCC constructed using two 
constituent punctured RSC codes. 

Puncturing of a rate- 1/2 RSC to obtain a higher rate RSC 
is represented by an 2 x M matrix as follows: 



(1) 



where M is the puncturing period and pi,m G {0, 1}, with 
i — and m—1,.... M. For „ — 0, the corresponding 
output bit is punctured. The puncturing pattern P for the 
rate- 1/2 encoder consists of the puncturing vector for the 
systematic output sequence and the puncturing vector V z for 
the parity-check output sequence. 

It was shown in [2] and [3] that performance bounds for 
a PCCC can be obtained from the transfer functions, or 
equivalently the weight enumerating functions (WEFs), of the 
terminated constituent RSC codes. A WEE provides all paths 
of length N that start from the zero state, can remerge with 
and diverge from the zero state more than once, and terminate 
at the zero state. 

More specifically, the conditional WEE (CWEF) of a 
punctured convolutional code C, denoted as A'^ {w, U, Z), 
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assumes the form [3] 



(2) 



where „ ^ is the number of codeword sequences composed 
of a systematic and a parity-check sequence having weights u 
and z, respectively, which were generated by input sequences 
of a given weight w. The overall weight of a codeword 
sequence is u + z. 

The input-redundancy WEF (IRWEF), A^'{W,U,Z), 
provides all codeword sequences for all possible values of 
input information weight, and is related to the CWEF as 
follows [3] 

A'^'(Vy,f/,Z) = ^yl^'(iu,(7,Z)W'". (3) 

w 

A relationship between the CWEF of a PCCC and the 
CWEFs of the constituent codes, €[ and C2 respectively, 
can be easily derived only if we assume the use of a 
uniform interleaver of size N, an abstract probabilistic concept 
introduced in [3]. More specifically, if A'^'^{'w,U, Z) and 
A'^'2{w,U, Z) are the CWEFs of the constituent codes, the 
CWEF of the PCCC, A{w, U, Z), is equal to 

{w, U, Z) ■ A^'-^ (w, U = 1,Z) 



A{w, U, Z) = 



N 



(4) 



The systematic output sequence of the second constituent 
encoder is not transmitted, therefore it does not contribute to 
the overall weight of the turbo codeword sequences, so it is 
eliminated by setting U=\ in A'^''^{w,U, Z). The IRWEF of 
the PCCC, A{W, U, Z), can be computed from the CWEF 
A{w^ U, Z), in a manner identical to (O. 

The input-output weight enumerating function (lOWEF) 
provides the number of codewords generated by an input 
sequence of information weight w, whose overall weight is 
d, in contrast with the IRWEF, which distinguishes between 
the systematic and the parity-check weights. For the case of a 
punctured PCCC, the corresponding lOWEF assumes the form 



B{W,D) = 



(5) 



where the coefficients B 



w.d 



can be derived from the 



coefficients Aw^u,z of the IRWEF, based on the expression 



Bu 



= E 



(6) 



The lOWEF coefficients i?t„ ^ can be used to determine 
a tight upper bound, denoted as P^, on the BEP Pb, for 
maximum-likelihood (ML) soft decoding for the case of an 
AWGN channel, as follows [3] 



(7) 



where P{w) is the union bound of all error events with 
information weight w, and is defined as 
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Fig. 1. Contribution to the union bound for rate-1/2 S-PCCC(1, 17/15, 17/15) 
using an interleaver of size A'^ = 1,000 and Af = 10,000 bits 



where R is the rate of the punctured turbo code. 

In [4], Benedetto et al. investigated the performance of 
rate-1/3 PCCCs and observed that the union bound P(wrnin) 
of all error events with the lowest information weight Wmin^ 
becomes dominant as the interleaver size N increases. Owing 
to the structure of an RSC encoder, the minimum information 
weight of a terminated RSC code is always equal to two, i.e., 
Wmm = 2. Consequently, the overall performance bound P^ can 
be approximated by P{2), when a large interleaver size is used. 
The same trend is also observed in the case of punctured turbo 
codes. The contribution, as a percentage, of P{2) and P(3) to 
P^ is illustrated in FigH] As an example, rate-1/2 S-PCCC(1, 
17/15, 17/15) is considered, using a uniform interleaver of 
size either iV= 1, 000 or iV = 10, 000. It is apparent that P{2) 
becomes the dominant contribution over a broad range of BEP 
values, as the interleaver size increases. 

We see from (|8]l that P{2) depends heavily on the minimum 

weight of the turbo codeword sequences, commonly known 

c' 

as free effective distance cifree.eff [4]. We use the notation d^^^ 

to denote the minimum weight of the codeword sequences 

c' 

generated by the first constituent encoder and zj'^^^ to denote 
the minimum weight of the parity-check output sequences, 
generated by the second constituent encoder. In both cases an 
input sequence of information weight 2 is assumed. Therefore, 
the free effective distance of a PCCC can be expressed as 



^free.eff 



(9) 



The free effective distance is the most significant parameter 
that influences the PCCC performance. The constituent 



encoders should be chosen to maximise d. 
consequently dfi-ee.eff- 



and and 



III. Computing the Upper Bound Approximation 

In order to compute P{2) and thus obtain a good 
approximation to the overall performance bound P^, we only 
need to calculate the CWEF of each constituent code for w = 2, 
i.e., A'^'^ (2, U, Z) and A^'^ (2, U, Z). 

Both CWEFs could be obtained by brute-force, i.e., input 
all possible sequences of weight 2 to each constituent encoder 



and group the output codeword sequences according to their 
systematic and parity-check weights. Although this approach 
is conceptually simple, it is extremely time-consuming, 
especially when a large interleaver size is used. 

The techniques proposed in [7] and [10] are more complex 
but less time-consuming. They both use the state diagram of 
a parent RSC code and introduce the puncturing patterns to 
obtain the full CWEF of the corresponding punctured RSC 
code. However, for large interleaver sizes and puncturing 
patterns with a long period, complexity becomes a prohibiting 
factor for the implementation of either approach. 

In this section we use the properties of the trellis structure 
of RSC codes to express the CWEF, for w = 2, of an 
RSC encoder as a function of its memory size, generator 
polynomials, and puncturing pattern. Consequently, derivation 
of the state equations and computation of the full transfer 
function of each constituent code, required in [7] and [10], 
is not necessary. Hence, PCCCs using both a large interleaver 
and a long puncturing pattern can now be easily supported. 

A. Unpunctured Rate-1/2 RSC Encoders 

A rate-1/2 RSC encoder, C, is characterised by its 
feedback and feedforward polynomials, Gn{D) and Gf{D) 
respectively. The degree of each polynomial is equal to the 
memory size ly of the encoder. A hypothesis commonly made 
[1], [4] so as to facilitate analysis of RSC codes is that Gf{D) 
is a monic function and that the initial state of the encoder is 
the zero state, for every input sequence of length N. 

Input sequences of weight 2 force the trellis path to diverge 
from the zero state and re-merge with it, after a number of 
time-steps. More specifically, the input sequence will change 
the state from to 2"^^, when the first non-zero bit is input to 
the encoder, as it is illustrated in Figl2] For as long as a trail 
of zeros follows the first non-zero input bit, the RSC encoder 
behaves like a pseudo-random generator, with the same state 
transitions being repeated every L time-steps, where L is the 
period of the feedback polynomial. In order for the path to 
re-merge with the zero state, the second non-zero bit should be 
input to the encoder when state 1 is reached, i.e., after kL + 1 
time-steps, where fc = 1, 2, . . . , [(A^-l)/iJ and [(iV-l)/Lj 
is the integer part of {N—l)/L. Furthermore, as it is depicted 
in Fig|2] when a non-zero input bit causes the path to diverge 
from or re-merge to the zero state, both the systematic and the 
parity-check outputs give a logical 1 . Therefore, if z^g^.g is the 
parity-check weight due to the transitions of the encoder from 
state 2"^^^ to state 1, the overall weight z of a parity-check 
sequence can be expressed as 

z(fc) = fczi, + 2, forfc = l,2,...,L(iV-l)/Lj. (10) 

Note that the state sequence during the transitions from state 
2"^^ to state 1 and, consequently, the value of ^fore' depend on 
the selected feedback polynomial. The minimum parity-check 
weight can be derived from z{k) by setting k — l, i.e.. 



Based on © and the CWEF for u; = 2, A'^ (2, U, Z), of 
the rate-1/2 RSC code when no puncturing is applied, assumes 
the form 

L(JV-1)/LJ 

A^i2M,Z)^ J2 42,.wt/'^^''\ (12) 
fc=i 

where ^2 2z(fc) number of codeword sequences with 

parity-check weight z{k), given by 

42,.W ^N-kL. (13) 

When the feedback polynomial, Gu{D), of an RSC encoder 
is selected to be primitive, the encoder visits all possible 2"^ — 
1 states with a maximum period of i = 2"^ — 1 time-steps 
[11], if the information weight of the input sequence is 2. As 
pointed out in [12], maximization of L increases the length 
of the shortest weight 2 input sequence, therefore increasing 
the chance of achieving a high weight z^^^. and, consequently, 
-^min- exact cxpression for z^^^^ can be derived based on the 
properties of pseudo-random sequences [11] or the analysis in 
[4], i.e., 

^core — 2 , (14) 

provided that Gr{D) ^ Gf{D). 

Since zf^^^ only depends on the memory size of the encoder, 
so does the CWEF of each constituent code, A*^ {2,U, Z) 
and, consequently, the union bound of all error events with 
information weight 2, P(2). Therefore, the performance of a 
rate- 1/3 PCCC, using a large interleaver size, mainly depends 
on the memory size of each constituent RSC encoder and not 
the underlying code, provided that the feedforward polynomial 
of each RSC encoder is different from the feedback primitive 
polynomial. 

B. Punctured RSC Encoders 

Rates higher than 1/2 can be achieved using a 2 x M 
puncturing pattern P on a parent rate-1/2 RSC encoder C. At 
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Fig. 2. Trellis diagram for codeword sequences of information weight two. 
Dashed lines con'espond to paths generated when the rate-1/2 RSC encoder 
operates as a pseudo-random generator A pair Xi / yi next to a branch denotes 
the input and parity-check output bits, respectively, generated at the end of a 
state transition. 



a time step i (0 < i < N), the weights of the systematic and 
parity-check output bits of the punctured encoder C will be 
Xi ■pi,m and Ui ■p2.m, respectively, where Xi, yi are the output 
bits of the parent rate-1/2 encoder and pi^m, P2,m are the 
elements of column-m (1 < m < M) of the puncturing pattern 
P. Note that, owing to the systematic nature of the encoder, 
Xi also represents the input information bit. The relationship 
between m and i is 



rem(i + l,M), 



(15) 



where rem(i + 1, M) denotes the remainder from the division 
{i + 1)/M. Since the period of P is M, its elements are 
repeated in such a way that pi^™ = Pi,(m+jM) and p2,m = 
P2,(m+iM)' where j is a non-negative integer 

In order to compute the CWEF of the punctured RSC for 
information weig htw = 2, i.e., A'^'(2,C/,Z), we need to derive 
an expression for the weight of the systematic and parity-check 
output sequences. Although information sequences with w — 2 
generate paths of length kL + 1, we first consider paths of 
length L+1, i.e., fc = 1, for simplicity. The weight u{k — l, m) 
of a systematic sequence, whose path diverges from the zero 
state when pi^m is active, is given by 



w(fc=l,m) =Pl,.m +Pl,(m+L), 



(16) 



since the two non-zero bits occur at the very beginning and at 
the very end of the path. Similarly, the weight z(fc = 1, m) of 
the parity-check sequence, whose path diverges from the zero 
state when P2,m is active, assumes the form 



z(fc = l,m) =P2,m + ^c"re^ +P2,(m+L), 



(17) 



since the parity-check bits at the beginning and at the end 
of the path are non-zero, while the weight of the remaining 



path is 
calculate z{k 
to derive 



as it is illustrated in Fig|3] In order to 
l,m) for every value of m, we first need 
2:£,„, . . . , 2*^g, by applying the M circularly 



^core 7 core ' 



shifted versions of the puncturing vector [p2.i, • ■ • ,P2,j\/] to the 
corresponding output parity-check bits of the parent rate-1/2 
RSC encoder, i.e. 



L-l 

{Vi ■P2,(i+va-l) 

i=l 



(18) 



If we extend our analysis to codewords associated with paths 
of length kL+1, we obtain the generic expressions for u{k, m) 
and z{k,ni) as follows 



u{k, m) = pi^rn + Pl,{yn+kL) 
fc-1 

Z(fc,m) = P2,m + ^ 2c"ore^^^^ +P2,(m+fcL), 



(19) 
(20) 



where z^,^-'^ = z^^g, due to the periodicity of the puncturing 
pattern. Since any codeword sequence, generated by an input 
sequence of weight 2, can be described by a polynomial 
jju{k,m) ^z{k,m) f^j. ^ gjven k and m, the summation of all 
polynomials of the form jj^i^^,™.) ^z{k,m) ^^^^ ^jj possible 
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Fig. 3. Trellis diagram for the weight calculation of parity-check sequences 



Union bound P, W=100 

— e — Bound approximation P(2), W=100 
- - Union bound P, W-1 ,000 
-c^ Bound approximation P(2), W=1 ,000 

Union bound P, W=1 0,000 

- e - Bound approximation P(2), W=1 0,000 




Fig. 4. Exact union bounds and their approximation for a rate-1/2 
S-PCCC(1, 17/15, 17/15) using interleavers of various sizes 



values of k and m will give the CWEF, A^' (2, U, Z), of the 
punctured RSC code 



L(Af-l)/LJ M 

A'^'{2,U,Z)^ ^ ^f'„C/"('=^'")Z^('=''"\ (21) 



k—l m— 1 



where A'l is the total number of codeword sequences with 
systematic weight u{k,m) and parity-check weight z{k,m). 
Coefficients can be easily derived if we observe that 

there are N — kL codeword sequences of length kL + 1 each. 
The codeword sequences are grouped into AI groups, whose 
members share the same weights u{k,m) and z{k,m). Thus, 
the number of codeword sequences in the m-th group is given 
by 



k,m 



_ M _ 
N-kL 



if rem ((TV- fcL),M) <m 



+ 1, 



otherwise. 



(22) 



Using (|2TJ, we can accurately and efficiently derive P{2), 

1. e., the probability of all error events with information weight 

2, which is a good approximation of the union bound P^, for 
a large interleaver size. In the example shown in FigH] we see 
that P{2) closely matches P^", when the interleaver reaches 
the size of A^=10,000 bits. 



IV. Case Study: Pseudo-random Puncturing 

In this section we consider constituent RSC encoders 
employing primitive feedback polynomials, therefore the 
period L assumes the maximum value of 2"^ — 1. Furthermore, 
we assume that the elements of the puncturing vector Pz, for 
the parity-check output, form a pseudo-random sequence of 
period A/ — L, generated by the same primitive polynomial as 
that of the RSC encoder 

Since the puncturing period M is equal to the period L of 
the feedback polynomial, u{k, m) and z{k, m) are reduced to 

u{k, m) = u{m) = 2pi^m, (23) 

z(fc,m) = fcz,'^+i + 2p2.rn. (24) 

Calculation of z(fc, m) and, consequently, (2, U, Z), 
requires knowledge of the L values of z^^.^. However, the 
assumption of pseudo-random puncturing can further simplify 
the computation of z{k,ni). 

A. Derivation of the Minimum Weight Values 

In order to express in a more compact form, we first 
need to consider the autocorrelation function (f){j) of a polar 
sequence, which is defined as [13] 

L 

Hj) = - l)(22/.+j- - 1) (25) 

1=1 

where j/i = {0, 1} is the output of the parent rate- 1/2 RSC 
encoder at time-step i for an input sequence of information 
weight 2, and < j < L. The parity-check sequence 
generated during the time period from i = 1 until i = L is 
a pseudo-random sequence, provided that the encoder does 
not return to the zero state. In this case, the autocorrelation 
function can be reduced to [11], [14] 

2'' - 1, if j = 
-1, ifl<j<L. 



m = 

Combining dZSl l and ( |26] | we find that 



(26) 
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if j = 

if 1 < j < L. 



(27) 



Since the puncturing vector for the parity-check bits, ~ 
[P2.i\^ is also a pseudo-random sequence generated by the same 
primitive polynomial Gr{D), such that P2,i+i = Vi, we can 
rewrite expression dZTl i as follows 



L 

E 

i=l 



{Vt ■P2,{i+m)) 



2^-2 



if TO = 1 

if 2 < TO < i 



(28) 



where j was replaced by m— 1. Due to the structure of the RSC 
encoder, the last bit of the parity-check sequence is always 
zero, i.e., i/l^O, therefore 

if m = 1 
if 2 < m < L 



and that L = 2"^ — 1, the weight of a parity-check sequence 
assumes the form 



z(fc, to) 



kT-^, 
fc2''-2 



2p2 



if TO = 1 

if 2 < TO < 2"- 1 



(30) 



where 2"^^^ elements of the puncturing vector P^ are equal to 
1 and the remaining 2"^^—! are equal to 0, since the elements 
of Pz form a pseudo-random sequence [11], [14]. 

The minimum weight of the parity-check sequences, 
can be expressed as 

zj;;„ = min {z(fc = 1,to)} = i ^^'^2 f"! !^ ^ o (31) 



7n—l...L 



for > 2 



whereas the minimum weight of the codeword sequences, d!^^ 
assumes the form 

d^i,, = min {u{m) + z{k — 1,to)} . (32) 

va—l...L 

As in the case of rate- 1/3 PCCCs, we conclude that when 
a large interleaver is used, the performance of a PCCC whose 
parity-check sequences were punctured using pseudo-random 
patterns, mainly depends on the memory size of the constituent 
RSC encoders, and not the exact underlying codes. 

B. Example Configurations for Rate-1/2 PCCCs 

In order to maximize the minimum weight of the codeword 

c' 

sequences, d^^^, generated by the first constituent RSC encoder 
of a PCCC, we can set the puncturing vector for the systematic 
output, P[/ = [pi,m], to be the complement of the puncturing 
vector for the parity-check output P^ = [P2.m], i-e-, Pi,m = 
P2.m- This configuration prevents u{m) and z{k,m) from 
assuming the smallest values at the same time. Therefore, 
expression {3% becomes 



"min 



-)iy-2 



(33) 



^iVr ■P2,(r+m)) = Cre^ = \ j-^-z' oV 1 ^ T (29) 

so z^jg is now a function of the memory v of the RSC encoder, first constituent encoder is given by (l33]l, while the minimum 
Having in mind that p2,L+i = since ul = 0, and so is P2,i, parity-check weight z^^?^ for the second constituent encoder is 



for v>2. 

A code rate of 1/2 can be achieved, if the parity-check 
output of the second RSC encoder is not punctured. In that 
case, z^jjj can be derived from (fTTI) and (fT4l l. The free effective 
distance of the corresponding PS-PCCC assumes the form 

4-ee,eif = 4 + 3(2''-2). (34) 

We refer to this example configuration as "Pseudo A". 

If our objective is to obtain a turbo code whose BEP 
performance quickly converges to the union bound but 
experiences a high error floor, we need to increase the number 
of transmitted systematic bits [10], [15], [16]. The parity-check 
output of both the first and the second constituent encoder is 
punctured using the same vector V z- Bearing in mind that 
'Pjj is taken to be the complement of V z, we need to replace 

all but one of the O's in Vu with I's, in order to achieve a 

c' 

code rate of 1/2. The minimum codeword weight d^^^^ for the 



TABLE I 

Puncturing patterns for rate-1/2 PCCC(1, 17/15, 17/15) 





Pseudo A 


Pseudo B 


Litt A 


Litt B 


Vector for Sys. Output 


[1000101] 


[1111101] 


[0010] 


[11] 


Vector for 1st Par.Output 


[0111010] 


[0111010] 


[1101] 


[10] 


Vector for 2nd Par Output 


[1111111] 


[0111010] 


[1111] 


[01] 



given by ( [3T1 ). The summation of the two minimum weights 
yields the free effective distance of the PS-PCCC 



^free,eff 



5, 
2H 



for v^2 
for V >2. 



(35) 



We refer to this example configuration as "Pseudo B". 

The particular puncturing patterns of each example 
configuration for the case of PCCC(1, 17/15, 17/15) are 
presented in Table U The configuration denoted as "Litt 
A" achieves a very low error floor and it was obtained 
through exhaustive search using [10], whereas "Litt B" is the 
conventional approach for obtaining rate-1/2 turbo codes. 

C. The Benefits of Pseudo-random Patterns 

Good punctured PCCCs can only be found by means of 
an exhaustive search among all possible patterns of a specific 
puncturing period M. The selection of a good pattern is not 
intuitive, since it can lead to catastrophic puncturing [15], 
i.e., rf^jij = 0, or semi-catastrophic puncturing, i.e., = 
for some values of to, of a constituent code C Furthermore, 
calculation of d^jj^ and z'^^^ requires prior knowledge of the 
M values of z™^^. 

The selection of a pseudo-random puncturing pattern 
guarantees that z™^^ > 0, and consequently, d^j^ > 0. Moreover, 
z™jg can be expressed as a function of the memory size v 
of C, permitting the immediate derivation of the minimum 
weights that characterize the PCCC. For the given puncturing 
rate of the parity-check output, the minimum value of 
is maximised, and so is z^j^, due to the properties of 
pseudo-random sequences. 

In Fig|5] we have plotted the performance of all four 
rate-1/2 PCCC(1, 17/15, 17/15) configurations, presented in 
Table U We observe that "Pseudo B" slightly outperforms the 
conventional "Litt B" configuration, while the performance of 
the PCCC based on the easy to derive "Pseudo A" pattern is 
close to the performance of the PCCC based on the "Litt A" 
pattern, obtained through exhaustive search. 

V. Conclusions 

We presented a simple approach to calculate the CWEF 
of punctured RSC codes, for input sequences of minimum 
information weight, which facilitates the approximation of 
the upper bound to the BEP, for punctured PCCCs using 
large interleaver sizes. Our technique offers the advantage of 
simplicity and reduced complexity, compared to time-hungry 
approaches, such as brute-force, or the more complex methods 
developed in [7], [10]. 

Furthermore, we considered pseudo-random puncturing 
patterns as a case study for our technique and we demonstrated 



Bound approximation P(2), Pseudo A 
Bound approximation P{2), Pseudo B 
Bound approximation P{2). Litt A 
Bound approximation P(2), Litt B 




Fig. 5. Bounds for various rate-1/2 PCCC(1, 17/15, 17/15) configurations, 
using an interleaver of size 1 ,000 bits 



that they prevent catastrophic or semi-catastrophic puncturing 
and facilitate the calculation of the minimum output weights 
of a turbo encoder, which characterize the performance of 
PCCCs. We concluded that pseudo-random puncturing could 
be used to obtain rate-1/2 PCCCs exhibiting low error floors, 
while specific puncturing patterns that achieve either a lower 
error floor or quicker convergence to the ML performance 
bound, could be determined by a subsequent search. 
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